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Abstract 

This paper presents a numerical study of the thermal performance of fins mounted on the bottom wall of a 
horizontal channel and cooled with either pure water or an Al 2 0 3 -water nanofluid. The bottom wall of the channel 
is heated at a constant temperature and cooled by mixed convection of laminar flow at a relatively low 
temperature. The results of the numerical simulation indicate that the heat transfer rate of fins is significantly 
affected by the Reynolds number (Re) and the thermal conductivity of the fins. The influence of the solid volume 
fraction on the increase of heat transfer is more noticeable at higher values of the Re. 
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Background 

Recent years have eyewitnessed a blossom flourishing in 
the evolvement of electronics, communications, and auto- 
computing industries, and this bearing is irrefutably con- 
tinuing in this century. The cooling of electrical, mechan- 
ical, and electronic components has become troublesome 
in today's fast-growing technologies. Inasmuch as the sig- 
nificance of heat exchangers in tremendous engineering 
applications, the subject of potential heat transfer enhance- 
ment in these devices has received sizeable attention in 
practice and research. On account of the fact that the 
consistency of the electronic components commodiously 
increases, conspicuous lack of heat transfer enhancement 
both in macro- and microscales of channels is realized. 
Encountering a fluid flow by utilizing transverse surfaces in 
a channel is a prevalent method that is used to intensify 
the rate of heat transfer from heated surfaces. 

Alamyane and Mohamad [1] studied the forced convec- 
tion heat transfer in a channel with extended surfaces. The 
effects of the Reynolds number (Re) and the fin height and 
spacing on the fluid flow and the heat transfer were exam- 
ined. Yang et al. [2] simulated the forced convection in a 
parallel plate channel. Constant temperature was consid- 
ered in both upper and lower walls, and a transverse object 
was located at the lower channel wall. The effects of the 
Reynolds number, the thermal conductivity ratio of the 
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fluid, and the fin profile area on the fluid flow and the heat 
transfer rate were analyzed. The study results showed that 
the heat transfer enhancement with an increment of the 
Reynolds number and the thermal conductivity ratio of the 
fluid at various fin profiles. Yang et al. [3] numerically inves- 
tigated the effect of mix convection heat transfer in an 
inclined parallel plate channel with a transverse object at 
the bottom wall. In this research, the effects of thermal con- 
ductivity, Reynolds number, the fin profile, and the channel 
inclination on the heat transfer rate at various Richardson 
numbers were examined. They discovered that the ace 
aspect ratio of the fin was related to the fin with utmost 
heat transfer at various Reynolds and Richardson numbers. 

Young and Vafai [4] observed the impact of controlling 
parameters on the cooling of heated channels with 
mounted objects. Concentrating on the effect of altering 
the dimensions of the object, the thermal conductivity, the 
heating method, and the Re was embraced. They deduced 
that the fluid flow and heat transfer are affected by the 
geometry and material of the object, and a correlation for 
the average Nusselt number was proposed as a function of 
the controlling parameters. 

Meinders and Hanjalic [5] experimentally investigated 
the effect of the cubes' arrangement on the turbulent 
fluid flow. They comprehended that the flow stream was 
affected by the distance between the objects owing to 
the fact of augmenting the flow velocity. Moreover, 
amelioration in velocity distribution and heat transfer 
than the staggered distribution case was found for flow 
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over inline cubes. Yan et al [6] experimentally investigated 
the influence of short surface-mounted objects at the top of 
a flat plate on the heat transfer enhancement. Scrutinizing 
was done on the effect of varies cross sections, spacing and 
numbers of objects, and the Reynolds number. They 
perceived that the heat transfer was incremented when the 
height of the object is comparatively equal to half of the 
channel height. 

In an experimental investigation by Yuan et al. [7], the 
heat transfer and friction characteristics of a channel 
which were attached by winglets were examined. Heat 
transfer from the channel was achieved to be noticeably 
augmented by using winglets in comparison with conven- 
tional channels with rectangular transverse objects. For a 
high Reynolds number, the heat transfer was enhanced by 
a factor of 2.7 to 6 times of the smooth channel. 

Utilizing nanofluids for the purpose of enhancing the 
heat transfer in thermal systems is another alternative 
technique [8]. The thermal performance of different types 
of nanofluids has been the subject of many recent studies 
on forced, natural, and mixed convection problems. Several 
explorations have studied natural convection of nanofluids 
in cavities [9,10]. They argued that the addition of 
nanoparticles in the fluid indisputably increase the natural 
convection heat transfer. 

Chein and Huang [11] analyzed the cooling of two silicon 
microchannel heat sinks with a water-Cu nanofluid. The 
heat transfer and fraction coefficients were based on the 
theoretical models and the experimental correlations. They 
realized that the heat transfer performance of microchannels 
was greatly improved when nanofluids were added into base 
fluid as coolants without any extra pressure drop. 

Recently, Santra et al. [12] numerically investigated the 
effect of water-Cu nanofluid through parallel plate channel 
in laminar forced convection. A cold nanofluid was sent 
through the channel, and the walls of the channel were 
isothermally heated. The effects of the Reynolds number 
and the solid volume fraction on the heat transfer were 
studied by considering the fluid to be Newtonian and non- 
Newtonian. They observed that the rate of heat transfer 
increased with an increase of the Reynolds number and the 
solid volume fraction. The increase in the heat transfer was 
approximately the same for both scenarios. 



The lattice Boltzmann method (LBM) is another numer- 
ical method that is often used to simulate flow problems. 
LBM have been used for more than 2 decades as an alter- 
native numerical technique. In LBM, it is intended to 
model fluids as a collection of particles, which successively 
undergo collision and propagation over a discrete lattice 
mesh. Several lattice Boltzmann models have been 
proposed for the incompressible Navier- Stokes equations. 
A collision model was proposed by Bhatnagar et al. [13] to 
simplify the analysis of the lattice Boltzmann equation, 
which leads to the so-called lattice BGK model. Remarkable 
efforts have been conducted by many researchers that made 
this numerical method more attractive for fluid dynamics 
modeling, e.g., [14,15]. For more details about LBM 
and its application, kindly refer to the aforementioned 
publications. 

Most of the researches cited above considered the heat 
transfer enhancement by adding either the fin or using 
nanofluids. The main objective of this study is to examine 
both of these effects on the heat transfer performance. In 
general, previous works were performed to investigate dif- 
ferent cases of nanofluid flow and heat transfer in channels 
with mounted objects by focusing on changing geometries, 
arrangement, and dimensions of the objects. However, 
more efforts are needed in order to optimize the controlling 
parameters for best heat transfer enhancement. 

Methods 

Problem definition 

The geometry of the problem is shown in Figure 1. A cold 
mixture of base fluid (water) and the nanoparticles (alu- 
mina) is forced to flow into a channel that is heated from its 
bottom and kept at a constant high temperature, while the 
top wall is insulated. The channel aspect ratio is fixed at L/ 
H = 15. The Prandtl number is taken as 7.02, and the Reyn- 
olds numbers are 10, 50, and 100, whereas the extended 
surfaces' height to space ratio 1/S is 0.2, and the ratio be- 
tween the objects' height to the channel's height 1/H is 0.2. 

The flow is assumed as Newtonian, laminar, two- 
dimensional, and incompressible. In addition, it is assumed 
that the cold mixture of base fluid (water) and the solid 
spherical nanoparticles (alumina) is in thermal equilibrium, 
and it flows at the same velocity as a homogenous mixture. 



Flow Direction 



1 



Channel length (L) 



Figure 1 A schematic plot of flow in a channel. 
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Numerical simulation 

The D2Q9 LBM model is used to simulate fluid flow in 
two-dimensional channel with uniform grid size of Sx x Sy. 
The lattice Boltzmann equation (known as LBGK equation) 
with single relaxation time can be expressed as [13] 



1 



fi(x + e/A£, t + At) -fi(x, t) 

which can be reformulated as 

fi(x + e/Af , t + At) = Off* 9 + [1 - (Dflfi 



(2) 



where o)f = ^ and as the single relaxation time of the 

fluid, fi represents the particle distribution function, is the 
particle streaming velocity, and f* q is the local equilibrium 
distribution function. For D2Q9 model is given by [8] 



l + Sei-u + ^.u) 2 -^ 2 



(3) 



where p is the density of the fluid and cd t is the weight 
function, which has the values of o)q = |, o)i — \ for t = 1 
to 4, and cot = ^ for i = 5 to 8. The macroscopic fluid flow 
velocity in lattice units is represented by u. In the LBM, 
the fluid macroscopic quantities such as density, p, and 
flow momentum, pu are calculated using the distribution 

8 8 

function yj, and given by p = ^^fi and pu = ^J&ifi, res- 

i=0 i=0 

pectively. The streaming speed for particles in coordinate 
(x and y) directions (i.e., 1 to 4, see Figure 2) can be 
expressed as = cos(tt/2 (i - 1)), sin(7r/2 (i - 1)), whereas 
particles in diagonal directions (i.e., 5 to 8 in Figure 2) have 
velocities of e t = V2(cos(tt / '4(2/ - 9)) , sin{n / '4(2/ - 9))); 
however, the particle in the lattice center is at rest and has 
no streaming speed, i.e., e 0 = 0. 

The thermal part is simulated using another distribu- 
tion function for the temperature. For instance, g is used 
to simulate the distribution function of the dependent 
variable (temperature) in the lattice Boltzmann equation, 
and an approach similar to that used to simulate the 
fluid flow is utilized to simulate the temperature 



distribution. In addition, the algorithm suggested by 
Succi [15] is adopted throughout this work. The kinetic 
equation for the temperature distribution function with 
single relaxation time is given by: 

gi(x + e/Af, t + At) - ft(x, t) = - ^ ( gi - g* q ) (4) 



ft 



(1) which can be written in the form 



gi(x + e,-Af, t + At) = co t g- q + [1 - G) t ] gl 



(5) 



Where g t represents the temperature distribution func- 
tion of the particles, g^ q is the local equilibrium distribu- 
tion function of the temperature, and o) t = ^, where f t is 
the single relaxation time of the temperature distribu- 
tion. Thus, the equilibrium distribution function of the 
thermal part is given by [15]: 



g eq = <p(x, t)cot 



1 + 



e, u 



C? 



(6) 



where, cp is the macroscopic temperature and Cj is the 
speed of sound. The diffusion coefficient can be obtained 
as a function of the relaxation time and given by ^ Q — |) 
. The macroscopic temperature is then computed from: 



(7) 



i=0 



A uniform lattice of 100 x 1,500 is used to perform all 
of the simulations. However, the number of lattices was 
doubled to test the grid dependency results. 

Since the inlet velocity of the flow is specified, the 
inward distribution functions should be computed at the 
boundary. In the D2Q9 model, the values of the distri- 
bution functions pointing out of the domain at the inlet 
boundary (i.e., f 3 , f 6 , f 7 in Figure 2) are known from the 
streaming step, and the only unknowns are (f lt f 5 , f 8 ) as 
well as the fluid density p. Following the work of Zou 
and He [16], the inlet density and the distribution 
functions can be obtained from: 

Pm = r^— \fo +/2 +/4 + 2(/s +fe +fi)\ (8) 

1 Ujyi 
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Figure 2 A schematic plot showing the thermal boundary conditions of the problem. 



Sidik et al. Nanoscale Research Letters 201 3, 8:1 78 
http://www.nanoscalereslett.eom/content/8/1 /1 78 



Page 4 of 8 



The unknown distribution functions are calculated 
using 



h =f3+^ pUx 



(9) 



An extrapolation scheme is used to simulate the outlet 
flow condition, which can be represented as fi(N xt t) = f t 
(N x - 1, t), i = 3, 6, 7. The bounce-back scheme is used 
to specify the boundary conditions on solid surfaces (no- 
slip boundary), in which the distribution functions 
pointing to the fluid are equal to those pointing out of 
the domain. The thermal boundary conditions for this 
case are given in Figure 2. For constant wall temperature 
(the lower wall temperature is constant), the unknown 
functions are obtained using the following equation [15]: 

gi = <Pwau(°>i + <*>i+*) - gi+2 for i = 2,5 and 6 (10) 

The left-hand boundary (channel inlet) is kept at a con- 
stant temperature (Dirichlet boundary condition) and set 
to a dimensionless value of zero. Therefore, the resulted 
equations of the unknown distribution functions on the left 
boundary are given by gi(0, t) = - g t + 2 (0, t), for i = 1, 5, 
and 8. For the adiabatic boundary condition, the gradient 
of the dependent variable normal to the boundary should 
be zero, i.e., d (p/d y = 0. The distribution functions are 
found to be in the following form [15]: 



(11) 



i=0 



i=0 



A second-order extrapolation similar to the one given 
in [17] is used to obtain the values of the unknown dis- 
tribution functions for the right-hand side boundary 
(channel outlet) as follows: 

gi(N x , t) = 2gi{N x -1,0- gi(N x - 2, f )for i 

= 3, 6, and 7. (12) 

The local Nusselt number (Nu x ) is computed using the 
following equation: 



Nu x 



L c (d(p/dy) y=Q 

Vwall ~ Wm 



(13) 



where L c is the characteristic length and (p wa ii is the wall 
constant temperature. The mean temperature cp m is 
given by: 



I 



ucpdA 



J udA 



(14) 



The effective density of the nanofluid is 

Pnf = (1 - <P)pf + <PPs 



(15) 



where cp is the solid volume fraction. The effective 
dynamic viscosity of the nanofluid given by Brinkman 
[18] is 




0 0 -F " i ~"' — i — i — i — i — i — i — i — i — i — i — i 
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Figure 3 Velocity and profiles at different cross sections. 
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Figure 4 Temperature profiles at different cross sections. 
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The thermal diffusivity of the nanofluid is 

OC n f — , v 

\P C p)nf 

The heat capacitance of the nanofluid is 

{P c p)nf = i 1 ~ <P)(p C p)f + (P C p)s 



(16) 



(17) 



k e ff is the effective thermal conductivity of the nanofluid 
and is determined using the model proposed by Patel 
et al [19]. For the two-component entity of spherical 
particle suspension, the model gives: 



Keff = Kf 



1 + + ck s Pe 



k f A f 



k f A f. 



(19) 



where k s and kf are the thermal conductivities of 
dispersed A1 2 0 3 nanoparticles and pure water. 



Pe 



Ctf 



(20) 



where u s is the Brownian motion velocity of the 
nanoparticles given by: 



2k b T 



(21) 



where kb = 1.3087xlO _23 JK _1 is the Boltzmann constant. 
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Figure 5 Velocity and temperature profiles at different Re. 
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Figure 6 Streamlines and isotherms for the Al 2 0 3 -water nanofluid and pure water at different Reynolds number. (A) Streamline plots at 
(a) Re = 10, (b) Re = 50, and (c) Re = 100. (B) Isotherm plots at Re = 10 and (a) cp = 0.0 and (b) cp = 0.05. (C) Isotherm plots at Re = 50 and (a) 
cp = 0.0 and (b) cp = 0.05. (D) Isotherm plots at Re = 100 and (a) cp = 0.0 and (b) cp = 0.05. 



Results and discussion 

Code validation and computational results 

For the purpose to ensure that the obtained results are 
proper and that the code is free of errors, a flow of cold 
air in a two-dimensional heated channel was taken as a 
benchmark test. Both upper and lower walls were 
heated. The comparisons were carried up between the 
dimensionless velocity and temperature fields at different 



locations in the channel as shown in Figures 3 and 4. 
The obtained results were found to be identical to the 
results of [20]. 

Figure 5 shows the effect of Reynolds on the 
temperature profiles at the same cross sections for 
Re = 10, 50, and 100. The figures depicted that the 
temperature profiles are less sensitive to the change in 
Reynolds compared to the velocity profiles. 
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Figure 7 Average Nusselt number for various Re. 



The effects of the Reynolds number and the solid 
volume fraction on the heat transfer, isotherms, and 
streamlines are studied. Figure 6 presents the stream- 
lines and the isotherms for the Al 2 0 3 -water nanofluid 
((p = 0.05) and pure water at different Reynolds number 
(Re = 10, 50, and 100). 

The streamlines show that as the Reynolds number 
increases, the vortices that are formed behind the fins 
become larger and stronger. This can be more clearly 
illustrated in Figure 5 where the horizontal velocity in the 
middle section between fins is presented. At Re = 10, the 
velocity is consistently positive. However, as the Reynolds 
number increases, the flow velocity becomes negative. 
This is an indication of flow reversal. The strong vortex at 
high numbers enhances the heat transfer from left 
face objects to right face objects and the wall between 
the two fins. This difference, however, becomes noticeable 
at higher Re. 

At low Reynolds numbers, the conduction is the 
dominating mechanism of heat transfer. Therefore, 
the isotherms stretch above the fins and take a large 
area in the channel. As Re increases, the convection 
becomes the dominating mechanism, and the strong 
cold inlet flow pushes the isotherms near the bottom 
wall. The comparison between the isotherms of the 
nanofluid and pure water shows that in each point of 
the channel, the nanofluid temperature is higher than 
the pure water. It is due to the nanofluid's higher 
thermal conductivity. 

The current investigation is wrapped with the analysis 
of the effect of the Reynolds number and percentage of 
nanoparticle volume fraction on the heat transfer 
enhancement in the channel. Figure 7 and Table 1 
display values of average Nusselt number at various 



Reynolds numbers and solid volume fraction from 0% 
to 5%. These figures demonstrate that the Nusselt 
number increases with the Reynolds number for 
values of volume fraction tested in the present study. 
For example, at Re = 100, in the addition of volume 
fraction of 5%, the average Nusselt number increases 
about 17%. High Reynolds number results in high 
energy transport through the fluid and cause irregular 
motion of nanoparticle. The higher solid volume 
fraction further stimulates the flow and contributes to 
higher Nusselt number as shown in the figure. The 
presence of nanoparticles also increases the rate of 
heat transfer by conduction mode through the flow. 

Conclusions 

LBM was applied to simulate forced convection heat 
transfer in two-dimensional channel including extended 
surfaces to investigate the effect of changing different 
parameters such as Reynolds number (10, 50, and 100) 
and nanofluid (A1 2 0 3 ) volume fractions (0.0, 0.03, and 
0.05). The results showed that as the Reynolds number 
increases, the rate of heat transfer also increases. The 
formation of vortices both in front and behind the 
objects enhances the heat transfer process. As the solid 
volume fraction increases, the heat transfer is enhanced 

Table 1 Average Nusselt number for various Reynolds 
number and solid volume fraction 



Reynolds 
number 


Average Nusselt 
number 


(p= 0.0 


<p = 0.03 


(p = 0.05 


Re = 10 


Nu ave 


2.712 


2.826 


2.965 


Re = 50 


Nu ave 


5.294 


5.683 


5.919 


Re = 100 


NUave 


10.252 


10.797 


11.109 
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for all values of the Reynolds numbers. This enhance- 
ment is more significant at high Reynolds numbers. The 
heat transfer rate of the fins increases with the thermal 
conductivity ratio of the fin to pure water. This enhance- 
ment has a finite limit. At this limit, the temperature at 
all surfaces of the fins approach the wall temperature. In 
this condition, the fins behave like constant temperature 
of heat sources. 
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